Stability of local quantum dissipative systems 



Toby S. Cubitt* 1 ' 2 , Angelo Lucia^ 1 , Spyridon Michalakis* 3 , and 
David Perez- Garcia^ 1 

1 Departamento de Andlisis Matemdtico, Universidad Complutense 
de Madrid, Plaza de Ciencias 3, Ciudad Universitaria, 28040 

Madrid, Spain 

2 Centre for Quantum Information and Foundations, DAMTP, 
University of Cambridge, Centre for Mathematical Sciences, 
Wilberforce Road, Cambridge CBS OWA, United Kingdom 
3 Institute for Quantum Information, Caltech, Pasadena, CA 91125, 

USA 



Abstract 

Open quantum systems weakly coupled to the environment are modelled 
by completely positive, trace preserving semigroups of linear maps. The 
generators of such evolutions are called Lindbladians. For practical and 
theoretical reasons, it is crucial to estimate the impact that noise or errors 
in the generating Lindbladian can have on the evolution. In the setting 
of quantum many-body systems on a lattice it is natural to consider 
local or exponentially decaying interactions. We show that in this case 
local observables and correlation functions are stable if the Lindbladian 
is frustration free, translational invariant (uniformly), has a unique fix 
point (with no restriction on its rank) and has a mixing time which scales 
logarithmically with the system size. These conditions can be relaxed to the 
non-translational invariant case, at the cost of requiring Local Topological 
Quantum Order. As a main example we prove that classical Glauber 
dynamics is stable under local perturbations, including perturbations in 
the transition rates which do not preserve detailed balance. 



Contents 

1 Background and previous work 2 

2 Stability of open quantum systems 3 

3 Setup and notation 5 

* tsc25@cam.ac.uk 
t anlucia@ucm.es 
' spiros@caltech.edu 
§ dperezga@ucm.es 



1 



4 Main result 7 

4.1 Assumptions for stability 7 

4.2 Stability 9 

4.3 Local observables vs. global observables 10 

4.4 Do we need all the assumptions? 11 

5 Toolbox for the proof 13 

5.1 Lieb-Robinson bounds for Lindbladian evolution 13 

5.2 LTQO and frustration freeness 15 

5.3 Local rapid mixing 15 

6 Proof of main result 17 

6.1 Step la: from truncated to local rapid mixing 17 

6.2 Step lb: from global to local rapid mixing 20 

6.3 LTQO and unique fixed point 22 

6.4 Step 2: from local rapid mixing to stability 24 

7 Glauber dynamics 26 

8 Conclusions and open questions 32 
Appendices 32 
Appendix A The non-stable example 32 



1 Background and previous work 

The physical properties of a closed many-body quantum system are encoded 
in its Hamiltonian. Theoretical models of such systems typically assume some 
form of local structure, whereby the Hamiltonian decomposes into a sum over 
interactions between subsets of nearby particles. Similarly, the behaviour of an 
open many-body quantum system is encoded in its Liouvillian. Again, this is 
typically assumed to have a local structure, decomposing into a sum over local 
Liouvillians acting on subsets of nearby particles. 

Crucial to justifying such theoretical models is the question of whether their 
physical properties are stable under small perturbations to the local interactions. 
If the physical properties of a many-body Hamiltonian or Liouvillian depend 
sensitively on the precise mathematical form of those local terms, then it is 
difficult to conclude anything about real physical systems, whose interactions 
will always deviate somewhat from theory. 

Quantum information theory has motivated another perspective on many- 
body Hamiltonians. Rather than studying models of naturally occurring systems, 
it studies how many-body systems could be engineered to produce desirable 
behaviour, such as long-term storage of information in quantum memories [8, 
15, 16, 38], processing of quantum information for quantum computing [7, 9, 
10, 29, 41], or simulation of other quantum systems which are computationally 
intractable by classical means [2, 4, 5, 22, 24]. Again, stability of these systems 
under local perturbations is crucial, otherwise even tiny imperfections will 
destroy the desired properties. This has been studied particularly in the context 
of self-correcting topological quantum memories, where one requires stability 
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also against local sources of dissipative noise, and the relevant quantity is the 
time it takes to introduce logical errors in the system. It has been known since 
[1, 8] that a self-correcting quantum memory with local interactions is possible 
in four spatial dimensions. With the breakthrough given by the Haah code [15], 
it now seems plausible to find such a self-correcting quantum memory in 3D. 

Recently, and partly motivated by the dissipative nature of noise, this ap- 
proach has been extended to open quantum systems and many-body Liouvillians. 
First theoretical [30, 47] and then experimental [3, 31] work showed that creating 
many-body quantum states as fixed points of engineered dissipative Markovian 
evolutions could be more robust against undesirable errors and maintain coher- 
ence for longer. Intuitively, there is an inherent robustness in such models: the 
target state is independent of the initial state. If the dissipation is engineered 
perfectly, the system will always be driven back towards the desired state. This 
idea can be used to engineer dissipative systems both for storing quantum infor- 
mation and for carrying out computation via dissipative dynamics. However, it 
does not guarantee stability against errors in the engineered Liouvillian itself. 
Once again, stability against local perturbations - this time for many-body 
Liouvillians rather than Hamiltonians - is of crucial importance. 

In the case of closed systems governed by Hamiltonians, major recent break- 
throughs have given rigorous mathematical justification for believing in stability 
of the physical properties of many-body Hamiltonians. Starting with [6], it 
culminated in the work of [39] which proved that, under a set of mathematically 
well-defined and physically reasonable conditions, gapped many-body Hamiltoni- 
ans are stable under perturbations to the local interactions. 1 More precisely, in 
the presence of frustration-freeness, local topological quantum order, and local 
gap, the spectral gap of a Hamiltonian with (quasi) local interactions is stable 
against small (quasi) local perturbations (see Section 4 for a formal definition of 
these conditions) . The bound on the amount of imperfection tolerated by the 
system comes in terms of the local gap, the decay of the local topological order, 
and the strength (and decay rate) of the interactions. Furthermore, except for 
frustration freeness, which is a technical condition required in the proof, these 
conditions are in a sense tight. There exist simple counterexamples to stability 
if any one of the conditions is lifted. 

2 Stability of open quantum systems 

In this work, we study stability of many-body Liouvillians. Our main result 
shows that under analogous assumptions to those in [39], one can prove stability 
also in dissipative systems. 

However, although the result is analogous, the proof and even the definition 
of stability in the case of Liouvillians necessarily differ substantially from the 
Hamiltonian case. For Hamiltonians, the relevant issue is stability of the spectral 
gap. Via the quasi-adiabatic technique [18, 20], this in turn implies a smooth 
transition between the initial and perturbed ground states, showing that both 
are within the same phase. Note that the existence of a smooth transition does 

"'^Note that, in stark contrast to traditional perturbation theory, the perturbations considered 
here simultaneously change all the local interactions by a small amount. The perturbations 
are therefore unbounded, and standard perturbation theory is useless. It is the local structure 
of the Hamiltonian and the perturbations that ensures stability. 
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not imply that both ground states are close in norm, as the trivial example 
H = <8>|0)(0| vs. H = <g>(|0) + e|l))((0| + £<1|) shows. It does however imply a 
well-behaved perturbation in the expectation of local observables - such as order 
parameters - and correlation functions, which in most experimental situations 
are the only measurable quantities. 

For Liouvillians, we are interested in a definition of stability more related 
to the evolution itself, which accounts at the same time for both the speed of 
convergence and the properties of the fixed point. Here, we consider the strongest 
definition of stability: we want our systems (initial and perturbed) to evolve 
similarly for all times and all possible initial states. Thus, not only should the 
speed of convergence to the fixed points be similar, the fixed points themselves 
should be close and so should the approach to the fixed points. 

This definition is significantly stronger than stability of the spectral gap 
alone 2 , and is more directly relevant to the applications discussed above. As 
in the Hamiltonian case, the analogous trivial example shows that one cannot 
expect to attain such stability if we consider global measurements on the whole 
system. We therefore restrict our attention to local observables and few-body 
correlation functions. 

Furthermore, there are subtleties involved in extending this stronger definition 
of stability to dynamics with multiple fixed points. Rather than getting side- 
tracked by these, we defer consideration of multiple fixed points to a future 
paper, and restrict our attention to dissipative dynamics with unique fixed 
points. It is important to note, however, that we do not make any assumption 
on the form of the unique fixed point. In particular, we do not assume that 
it is full-rank (primitivity); our results apply equally well to Liouvillians with 
pure fixed points. (Pure-state fixed points are particularly relevant to quantum 
information applications, such as dissipative state engineering and computation.) 

A key technical ingredient in the stability proof for Hamiltonians is the 
quasi-adiabatic evolution technique [18, 20], which directly uses the fact that 
Hamiltonian dynamics is reversible. This is of course no longer true for Liou- 
villians, so we must use a different proof approach. We make use of the fact 
that evolution under a Liouvillian converges to a steady-state, together with 
dissipative generalizations [40] of the Lieb-Robinson bounds that are the other 
crucial ingredient in [39]. 

Among the systems which verify our hypotheses, one finds classical Glauber 
dynamics [37]. This immediately shows that Glauber dynamics is stable against 
errors. To the best of our knowledge, this is new even to the classical litera- 
ture. Given the importance of Glauber dynamics to sampling from the thermal 
distributions of classical spin systems [33, 37], we expect our results to have 
applications also to classical statistical mechanics. 

The paper is structured as follows. After setting up notation and basic 
definitions in the next section, we state our main stability result in section 4 
and discuss the assumptions it requires. In section 5 we prove various technical 
results used in the main proof, which is given in section 6. We apply these results 
in section 7 to the important example of classical Glauber dynamics, before 
concluding with a discussion of the results and related open questions in section 
8. 

2 Also note that by the recent work [44] it is not clear whether the spectral gap in Liouvillians 
is the relevant quantity for convergence questions. 
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3 Setup and notation 



We will consider a quantum system defined on a square lattice 3 A = [0,L- 1] D c 
TP or A = c T D . The ball centered at x e A of radius r will be denoted 

by b x (r). We will denote by N the number of sites in the lattice, i.e. N = L D . 
At each site x of the lattice we will associate one elementary quantum system 
with Hilbert space H x . Then for each subset T £ A the associated Hilbert space 
will be 

Hr = ®H x , 

and the algebra of observables supported on T is defined by 

Ar = <g)B(U x ). 

xer 

If Ti c T2, there is a natural inclusion of *4ri in Ar 2 by identifying it with 
.Ari ® 1- The support of an observable O 6 Aa is the minimal set T such that 
O = O' ® 1, for some O' € Ar, and will be denoted by suppO. We will denote by 
||-|| the Schatten p-norm over Aa- Where there is no risk of ambiguity, ||-|| will 
denote the usual operator norm (i.e. the Schatten oo-norm). 

A linear map T '■ Aa -*■ Aa will be called a superoperator. The support 
of a superoperator T is the minimal set r E A such that T = T' ® 1, where 
T' e B(Ar)- A superoperator is said to be Hermiticity preserving if it maps 
Hermitian operators to Hermitian operators. It is said to be positive if it maps 
positive operators (i.e. operators of the form 0*0) to positive operators. T is 
called completely positive if T ® 1 : Aa ® M n -*■ Aa ® M n is positive for all n ^ 1. 
Finally, we say that T is trace preserving if trT(p) = trp for all p € A a- For a 
general review on superoperators, see [49]. 

The dynamics of the system is generated by a superoperator C, which plays 
a similar role to the Hamiltonian in the non-dissipative case. The evolution will 
be given by the one parameter semigroup T t = e tc . The natural assumptions to 
make about T t are that it is a continuous semigroup of completely positive and 
trace preserving maps (CPTP, sometimes also called quantum channels). Such 
maps are always contractive, meaning that ||T(|| 1 ^. 1 cb ^ 1, where: 

IITII ll WT^KiX)^ 

ll T lli^i. C fc = su p|l T<g,1 nlli-i = su p su p — mi • (!) 

n n XeAA<SM„ II 111 

We will also be interested in the IHI^^^ cb norm of superoperators, which is 
defined as follows: 

IITII IIT^I II |T®l n (X)L 

\\ T \\ 00-^00,0, = SU P l T<g> 1 «L^=o = SU P SU P 



n XeAA®M n 
X*0 



11X11 



The relationship between ||- 1[ ^— cb anc ^ II"IIop->oo cb IS ^ ne following 

I-^l 1->1, cb = 11-^ Ilcx)^cx),c6 • 

We will denote IHIo^^ cb simply by ||-|| cb when there is no risk of confusing 
different completely-bounded norms. 

3 We restrict to square lattices for the sake of exposition. The results can be extended to 
more general settings, replacing the lattice [0, L — 1] with a graph equipped with a sufficiently 
regular metric. 
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Observation 3.1. As shown in [23], the supremum in equation (1) is reached when 
n is equal to the dimension of the space on which T is acting: if T : M n -*■ M n , 
then lT*l n \ M = 1^1^^. 

The generator C of such a semigroup is called a Lindbladian. All such 
generators can be written in the following general form, often called the Lindblad 
form [11, 34] (see [49]): 

Proposition 3.2. C generates a continuous semigroup of CPTP maps if and 
only if it can be written in the form: 

C(p) = i[p,H] + Y.hpL* ~ l{L*L JlP }, (2) 
j 1 

where H is an Hermitian matrix, {Lj}j a set of matrices called the Lindblad 
operators, [•,•] denotes the commutator and {•,•} the anticommutator. 

Since C is defined on a lattice A, it is natural to ask that it has some form of 
local structure. We will say that £ is a local Lindbladian if it can be written 
as a sum of terms, each of which is itself in Lindblad form and has bounded 
support: 

^ = SS ;C «( r )' supp£„(r) = b u (r). (3) 

u r 

Such a decomposition is obviously always trivially possible. We are interested 
in the cases in which the norms of C u (r) decay with r. Concretely, let us define 
the strength of interaction as the pair (J, /(r)) given by 

T II/. / M| P( \ ll^'«( 7 ')lll^l,cfc 

J = sup||£ u (r)|| 1 ^ lcfc , /(r) = sup . (4) 

u,r u <J 

We will restrict our attention to the case in which f(r) is compactly supported, 
which is usually called the finite range interaction case 

/(r) = Vr > k, for some finite k, 

or the case in which f(r) is decaying exponentially fast, 

f(r) = e~ Mr for some fi > 0. 

More generally, we will say that a function f(r) is fast-decaying if it goes to 
zero faster than any polynomial 4 , as r goes to infinity. In such case we say that 
C has quasi-local interactions. Quasi-local interactions verify a Lieb-Robinson 
bound (see section 5.1) with slower decay in space than the previous classes of 
interactions. For the sake of simplicity of the exposition, we will not go into such 
generalizations. Nonetheless, the proof presented in this paper applies, with the 
obvious modifications, to this setting. 

As shown in [48], from the spectral decomposition of C (and T t ) one can define 
two new CPTP maps which represent the infinite-time limit of the semigroup 
T t . We will denote by the projector on the subspace of stationary states 
(fixed points), and by the projector on the subspace of periodic states. They 



4 Or more generally, if it decays faster than a polynomial of some degree, which depends on 
the geometry of the lattice. 
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correspond, respectively, to the kernel of L and to the eigenspace of purely 
imaginary eigenvalues of C, which we denote Tc and Xc, respectively. Both the 
subspaces are invariant under T t : in particular, T t acts as the identity over Tc-, 
while it is a unitary operator over Xc . Note that both subspaces are spanned by 
positive operators (i.e. density matrices). 

As we are going to exploit the local structure of C, we will often be interested 
in the restriction of £ to a subset of the lattice. Given A c A, we define the 
truncated, or localized, generator: 

C A = Y. C u(r). (5) 

6„(r)c/t 

We will use superscripts to indicate truncated versions of the CPTP maps 
generated by Ca, so that T t , and indicate respectively the semigroup 
generated by the truncated generator, the projector on its fixed point subspace, 
and the projector on the periodic subspace of e* A . 

Definition 3.3. As usual, we define the spectral gap of T t to be the difference 
between 1, and the maximum modulus of the elements of the spectrum of T t 
strictly contained in the unit disk ( i. e. the spectrum of T t minus its peripheral 
spectrum). We say that T t is gapped if its spectral gap does not vanish with 
increasing system size. 

Definition 3.4. We say that C has a unique fixed point if Xc = = {Poo}- In 
other words, the spectrum of C has only one element with modulus 1; namely, 
the non- degenerate eigenvalue 1. 



4 Main result 

4.1 Assumptions for stability 

In order to prove stability of the generated dynamics, we will impose three 
conditions on the local Lindbladian. (These assumptions are analogues of the 
conditions for the stability of Hamiltonians imposed in [39] , and in some cases 
are direct generalizations of these to the dissipative setting.) 

Definition 4.1 (Frustation-Freeness) . We say that C is frustration free if for 
all Ac B c A, we have the following inclusion property for the periodic subspaces 
of the local Lindbladians: 

Xc b ^X Ca . (6) 

Topological quantum order, namely the property of some quantum states 
to be indistinguishable on a small scale, is a widely studied property of ground 
states of local Hamiltonians. In the dissipative setting there are no longer any 
ground states, but we can easily define the analogous concept for periodic states 
of Lindbladians. 

Definition 4.2 (Local Topological Quantum Order (LTQO)). 
Take a convex set A c A and let A(£) = {x e A| dist(a;, A) ^ £}. Given two states 
Pi 6 Xc A m> * = 1, 2 7 consider their reduced density matrices on A: 

Pi =^A(£)\APi, i = l,2. 
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We say that C has local topological quantum order (LTQO) if for each £ ^ 0: 

|K- ^1^(14) Ao09, (7) 

where Aq(£) is a fast- decaying function and k(\A\) is polynomial in \A\. 

In Hamiltonian systems, the spectral gap (the difference between the two 
lowest energy levels) plays a crucial role in a number of settings, from defining 
quantum phases and phase transitions [43] to understanding the entanglement 
and correlations present in the system [17, 19, 21] or analyzing its stability [6, 
39]. On the other hand, it is known that, for Lindbladians, the spectral gap alone 
is not sufficient to fully characterize the convergence properties of the dissipative 
evolution [27, 44]. Therefore, we will instead impose a more general requirement 
on the convergence of the dynamics. (It is not known - though it seems unlikely 
- whether this requirement depends only on simple spectral properties of £, i.e. 
properties depending on the eigenvalues - like the gap - and eigenvectors - like 
the condition number.) 

Definition 4.3 (Truncated rapid mixing). Define the contraction ofT t , 

V (T t )= sup dtr(Tt(p),TtT^p)), (8) 

trp=l 

where d% T denotes the trace distance d tr (p 7 o-) = | Hp- o~\ 1 . We say that a Lind- 
bladian C verifies truncated rapid mixing if for any x e A and any r ^ we have 
that 

r,(1*^)^k(r)e^, (9) 
for some 7 > and some k(r) > polynomial in r. 

Note that this definition implies that the global mixing time, i.e. the quantity 
t e = min{£ > | ?y(T t ) < e}, scales logarithmically with the system size. As 
we will see later, this seems to leave out interesting systems, for example the 
ones associated with dissipative state engineering of topological codes [47] . On 
the other hand a number of systems do verify this kind of condition, and new 
techniques are being developed that provide tools to show such convergence rates, 
like the use of quantum logarithmic Sobolev inequalities and hypercontractivity 
[25]. 

Up to now, we have made no assumption on how the local terms C u (r) at 
different lattice sites relate to each other, other than the very generic condition 
on the decay of their norms in r. Nor have we assumed anything about how C is 
defined for different system sizes. This very generic setting implies that equation 
(9) must be verified for all subsystems b u (r). For complex models this could 
be prohibitive. On the other hand, realistic models usually do not need such a 
generic formulation: the dynamics on different subsystems are typically related, 
or even identical in the most common case of full translational invariance. In 
section 6.2 we will give a precise definition of families of translationally-invariant 
local Lindbladians with arbitrary boundary conditions, which we call uniform, 

b (r) 

such that the truncated dynamics T t is similar to the global dynamics 
T t on a smaller lattice. For these Lindbladians, which in particular includes 
all translationally-invariant models with arbitrary boundary conditions, the 
assumption of definition 4.3 can be weakened. 
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Definition 4.4 (Global rapid mixing). We say that a Lindbladian C verifies 
global rapid mixing if 

ri(T t ) < fe(|A|)e-^; (10) 
for some 7 > and some k(r) > polynomial in r. 

This last assumption can be restated as a logarithmic scaling with system 
size of the mixing time t e for T t . Remarkably, if we assume that C has a unique 
fixed point for each system size (an assumption that will be required in the proof 
of the stability theorem) , then global rapid mixing is already sufficient to prove 
LTQO, for arbitrary boundary conditions. This is one of the first results we are 
aware of which gives conditions to prove LTQO. 

Theorem 4.5. Let C belong to a uniform family of local Lindbladians verifying 
global rapid mixing, frustration freeness, and having a unique stationary state. 
Then C verifies LTQO. 

The proof of this result, together with detailed definitions, will be given in 
sections 6.2 and 6.3. 

4.2 Stability 

With the required assumptions laid out, we can now state our main result. 

Theorem 4.6. Let C be a local Lindbladian with a unique fixed point, verifying 
frustration- freeness and either of the following conditions: 

(a) LTQO and truncated rapid mixing; or 

(b) C belongs to a uniform family verifying global rapid mixing. 

Let E = En E r -Eu( r ) be a sum of super- operators, which we will call the 
perturbation. Suppose the perturbation verifies the following assumptions: 

(i) \E^{r)\ ch ^ £f(r), where £ is a constant (the "strength" of the perturba- 
tion) and f(r) is a fast-decaying function. 

(11) £?*(r)[l] = 0. 

(Hi) e*' £+E ) is a contraction for each t^O. 
For an observable O a supported on Ac A, let 

O (t) = e tc '(O A ), and O^t) = e« c * +E *\0 A ). 
Then we have that for all t^ 

WOoM-o^w^cQadwoaWs, (11) 

for some c not depending on the system size, independent of t, and polynomial 
in \A\. 

Observation 4.7. The assumptions (ii)-(iii) on the perturbation E are satisfied if 
each C u (r) + E u (r) is a Lindbladian, but there are more general perturbations 
which are covered by the theorem. 
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Observation 4.8. Since we are free to choose an Oa with support on two non 
connected regions, we can apply theorem 6.11 to two-point correlation functions 
(or more generally fc-point correlation functions, for fixed k) and still obtain that 
the error introduced by the perturbation depends linearly on the strength of the 
perturbation (and not on its global norm). 

One of the tools already known in the setting of classical Markov chain 
[12-14, 37], and recently generalized to quantum dissipative systems [27], are the 
so called Logarithmic Sobolev inequalities (in short, Log Sobolev inequalities). 
Introduced in a different setting to study hypercontractivity of semigroups [28] , 
they provide the right asymptotic regime needed to verify the global rapid mixing 
hyptothesis: in fact, the existence of a system size independent Log Sobolev 
constant implies a logarithmic scaling of the mixing time, which is exactly what is 
required in Definition 4.4. Without going in the technical details of Log Sobolev 
inequalities, we summarize this fact in the following Corollary. 

Corollary 4.9. Let C belong to a uniform family of Lindbladians, verifying 
frustration-freeness and having a unique fixed point for each system size. If C 
verifies the Log Sobolev inequality with a system-size independent constant, then 
the system is stable, in the sense of theorem 4-6. 

A straightforward consequence of the techniques we use to prove stability, is 
that the fixed points local Lindbladians which satisfy the assumptions of theorem 
4.6 necessarily have fast decay of correlations. 5 

Theorem 4.10. Let C be a local Lindbladian satisfying either of the sets of 
conditions of theorem 4-6. Let Oa and Ob be operators supported on regions 
A and B, with \\Oa\\ = \\Ob\\ = 1, and call dAB be the distance between them. 
Define the correlation function 

cott(Oa,O b ) = \{O a O b ) - {Oa) (O b )\ , 

where (O) = tr(/9 00 0) is the expectation value of the operator O acting on the 
fixed point p x of C. Then we have that 

cott(O a , O b ) < P (\A\ + \B\)g{d AB ), 

where g(d) is a fast- decaying function and p(-) is growing poly nomially. Moreover, 
if C verifies global rapid mixing, or verifies LTQO with an exponentially decaying 
Aq, then g(d) is also exponentially decaying. 

4.3 Local observables vs. global observables 

The bound in equation (11) has a scaling on the size of the support of the 
observable B. Although the dependence is polynomial, for observables with 
large support the result is not useful. In most realistic experiments, only finite 
support observables and fc-site correlation functions can be measured, so this is 
not such an issue in practice. Nonetheless, one might ask more generally for a 
bound on snp p d tI (T rx ,(p),T^(p)) 1 where is the fixed-point projector for the 
perturbed C = C + E. 

5 As we were writing up this manuscript, a result similar to Theorem 4.10 was announced 
by Kastorayno et al. [26]. 
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However, this is not possible; the limitation to local observables is in some 
sense strict. There is no hope of finding such a bound for global observables, as 
the following simple example shows. 6 

Example 4.11. Consider N independent amplitude damping processes, with 
uniform rate 7 (which we can suppose w.l.o.g. equal to 1). This Lindbladian can 
be written as 

N 

£n = ^l-fc-i ® £1 ® lfc+l...JV, 
k=l 

where 

£ 1 (p) = |0)(l|p|l)(0|-i{p,|l)(l|} 

is an amplitude damping process on a single qubit, preparing the state |0). Then 
e tc N = ( e tC!j<$N hag gap jy 2 an d m i x i n g time of order O(logTV) [27]. The fixed 
point is of course the pure state |0 . . . 0)(0 . . . 0| . 

Now consider £f, a rotation of C\, which prepares the state \a$) = v 1 — e 2 |0) 

+e |1). We have - C\ || = O(e), but the new fixed point laoKaol 8 ^ ^ s a l m °st 
orthogonal to the original one, since the overlap between the two is 

(0 . . . 0|a . . . a ) = (0\a ) N = (1 - e 2 ) N/2 ~ e^' 2 -> as N ■* 00. 

This shows that in general the dependence on the support of the observable in 
equation (11) cannot be improved. 

4.4 Do we need all the assumptions? 

Are the assumptions of theorem 4.6 a reasonable set? We have just shown that 
we must necessarily consider local observables if we are to have meaningful 
bounds, but what about the other conditions? We will now present, in order, an 
example of each of the following: 

1. verifies frustration-freeness, truncated rapid mixing, global rapid mixing, 
and belongs to a uniform family, but does not have a unique fixed point 
and fails to verify LTQO; 

2. verifies frustration-freeness, LTQO, and has a unique fixed point, but fails 
to verify global rapid mixing, truncated rapid mixing, and does not belong 
to a uniform family; 

3. (presented in Appendix A) verifies frustration-freeness, LTQO, has a uni- 
que fixed point, and verifies global rapid mixing, but it neither verifies 
truncated rapid mixing, nor belongs to a uniform family. 

All these systems will be shown to be unstable. 

Example 4.12. Consider a ID chain composed of N 4- level systems, with an 
independent Lindbladian acting on each site, having the following Lindblad 
operators 

L 1 = |0)(1|, i 2 = |0)(3|, £ 3 = |2)(1|, i 4 = |2)(3|, 

^Indeed, all global stability results for quantum Linbladians we are aware of have a 
dependency on the total Hilbcrt space dimension [26, 45]. 
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and call 



4 1 

A)(P) = Y. L iP L i ~ ~{Pi L i L *}- 

i=l / 



The global Lindbladian Cm is given by applying Co independently on each site 
k = l...N; 

N 

£-N = Yj ® ^0 ® lfc+l....,Ar- 

fc=l 



Then we have that 



f|0)(0| + |2)(2|-2|i)(j| 




if » = J> ijj'e{l,3} 
if i = j, ije {0,2} 



,-[x{i,3}(«) + X{i,3}0')] if « * j- 



Diagonal states of the form \i)(i\ evolve according to the classical Markov 
process embedded in the Lindbladian, while off-diagonal elements \i){j\ evolve as 

T t (\i)(j\) = exp(-t[ X{ i,3}« + X{i,3}0')]) ■ 

This implies that T Co = span{|0)(0| , |2)(2| , |0)(2| , |2)(0|}. Since C has gap equal 
to 1, theorem 6.1 below implies that Cm verifies both truncated and global rapid 
mixing. Cm is frustration- free and forms a uniform family, but it does not verify 
either the unique fixed point condition, nor LTQO: clearly, the fixed points 
|0 . . . 0)(0 . . . 0| and |2 . . . 2)(2 . . . 2| can be distinguished locally. 
Consider now the following additional Lindbladian 



£oO) 



2 

N 



|0)(2|p|2)(0|--{p,|2)(2|} 



Then, we have: 



(A, + £b)(W0'l) 



|0)(0| + |2)(2|-2|i)(i| 

|(|o>(o|-K)(i|) 



-[x{i,a}(i) + X{i,3}U) + %^]\W\ 



if i = j, i,j 6 {1,3} 
if * = J = 2 
if i = j = 

if i * j- 



Again, this implies that Fc +£ = {|0)(0|}- Consequently Cm + £n verifies 
LTQO, frustration-freeness, and has a unique fixed point. It is not a uniform 
family, and it clearly does not verify truncated rapid mixing or global rapid 
mixing, as it is not even globally gapped: 



£ JV +£ A r(|2,0...0)(2,0...0|) 



(|0,0...0)(0,0...0|-|2,0...0)(2,0...0|). 



Analogously, Cm + £/v verifies the same conditions as Cm + £ni but the unique 
fixed point is now |2)(2|. 

All these three systems are unstable, since we can transform one into the other 
by applying a perturbation of order 0(l/iV), yet the fixed points of Cm + Sn 
and Cm +&m are locally orthogonal (while Cm has both of them as fixed points). 
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5 Toolbox for the proof 



The proof of theorem 4.6 will be done in two steps. First, we will introduce 
a new condition, called local rapid mixing, and prove that the hypotheses of 
theorem 4.6 imply local rapid mixing. The second step, independent of the first, 
will consist in proving that local rapid mixing and uniqueness of the fixed point 
imply the desired stability result. But first, we will need to introduce some useful 
tools. 

5.1 Lieb-Robinson bounds for Lindbladian evolution 

We first recall a generalization of Lieb-Robinson bounds to non-Hamiltonian 
evolution, due to [42] and [40], which we will use to derive a number of useful 
tools that allow us to approximate the support of an evolving observable with a 
finite set which grows linearly in time. 

Theorem 5.1. Let C = r C u (r) be a local Lindbladian satisfying the following 
condition for some positive fi, v 

E E I^tolli^i,c6Mr)|e^<oo VueA. (12) 

xeA r^dist(u,x) 

(The condition is verified automatically if C has finite-range or exponentially 
decaying interactions.) 

Let Ox be an observable supported on X c A, and denote by Ox(t) = T^(Ox) 
its evolution under C. Let K ; Ay -* Ay be a super- operator supported onYcA 
which vanishes on 1. Then it holds that 

\\K(0(t))\\ < ItfL^ IIOxl C(X,Y)e~^ x ^ (e* - l) , (13) 

where C(X,Y) = min(|X| , \Y\). v is called the Lieb-Robinson speed (or velocity) 
ofC. 

Note that if C has interaction strength (J, f(r)), then equation (12) can be 
replaced with 

E E /WIM0|e Mr <~<°° VueA. (14) 

xzA r2dist(u,x) J 

This implies that, if equation (14) is verified, different Lindbladians whose 
strengths can be uniformly bounded will have uniform Lieb-Robinson velocities. 
From now on, we will only consider Lindbladians which verify equation (12). 

Lemma 5.2 (Comparing different dynamics). Let C\ and C2 be two local 
Lindbladians, having strengths bounded by (J, f(r)). Let Y = supp(£i - £2) be 
the set on which the two differ. Consider Ox an operator supported on X c A, 
and call Oi(t) its evolution under Ci. Then it holds that 

WO^t) - 2 (t)\\ < ||A - C 2 \\ M>cb [[Ox! C{X,Y)e-^ x ^ 

(15) 
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Proof. Let h(t) = 0\{t) - 2 (t). Calculating its derivative, we obtain 

h'(t) = C[O x {t) - C* 2 2 (t) = C\h{t) + (£{ - C* 2 )0 2 {t). 

Since h(0) = 0, this differential equation for h(t) has solution 

t 

h(t) = 1 (t)-0 2 (t) = f e^^(C\-C* 2 )0 2 (s)ds, 
o 

giving us a useful integral representation for 0\(i) - 2 (t). From this, we obtain 
the estimate 

t 

\\Ht)H f \\(Cl-C;)0 2 (s)\\d S , 
o 

where we have used the fact that e fCl is a contraction with respect to Ij-jl^ for 
each t Z 0. 

We can now apply the Lieb-Robinson bound (equation (13)) to the previous 
estimate, to obtain 

t 

\\h(t)\\ < \\L X - C 2 \\^ cb \O x \ C(X,Y)e-^ x ^ f (e™ - 1) ds, 

o 

which implies the claimed bound. □ 

Lemma 5.3 (Localizing the evolution). Let Oa be an observable supported on 
finite Ac A. Denote by A (t) = T^(Oa) its evolution under a local Lindbladian 
C with strength (J,f(r)). Given r > 0, denote by O r (t) its evolution under the 
localized Lindbladian C A (r)> where A(r) = {x e A| dist(x,^4) ^ r}. 
Then it holds that 

\\0 A {t) - O r {t)\\ < \\O a \\ \A\ -(e vt - 1 - vt)g(r) (16) 

v 

where g(r) is an exponentially decaying function. 

Proof. First, let us decompose Oa(£) _ O r (t) as a telescopic sum 

o A (t)-o r (t) = Y,Oi + i(t)-o l (t). 

I2r 

Each element in the sum is the difference between the evolution generated by 
C-A{i) an d Ca(i+i)- Since they are restrictions on different subsets of the same 
global Lindbladian, it is easy to explicitly write the difference 

z+i 

£ a(i+i) - £- a(i) = E C u (i + i-6). 

<5=Odist(«,A)=c5 

We can then apply our previous integral repesentration for Oi + \(i) - Oi(t), and 
obtain: 

t 

A (t) - O r (t) = £ f e ^ c Mi^ (£* } - C A{l) )Oi(s) ds 

= E E E / e^M^c^l + 1 - S ) 0l (s) ds. 

V»rd=Q dist(M,A)=5 q 
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Now take norms, recalling that each e tCA <- 1 ) is a contraction: 

\OA(t)-o r (m<ZT, E / \\c:(i + i- 6)0^)^3. 

Vzr 5=0 dist(«,A)=5 g 

Since d = dist(A, b u (l+l-S)) = max{0, 26 -1-1} , we can apply the Lieb- Robinson 
bound (equation (13)) to each of the norms inside the integrals, and obtain 

f \\C u (l + 1 - a)0,(«) || ds < Jf(l + 1-6)\\O a || \A\ e-» d - t\ . 

We can group the terms in the sum by their distance from A. Denoting q(l) = 
\{u : dist(w, A) = l}\, we have 

\\0 A (t)-O r (t)\\ 



„vt _ _ 1 — ( + 1 



^\\0 A \\\A\j e v L^Y l Z9(S)f(l+l-S) + E q(6)f(l + l-5)e~^ 2S - l ~V 

V !>r6=0 5=l±i 

u 2 

p vt _ 01+ — 1 ' +1 

S IOaI |^1| J E E + 1 - 6)e^ l+1 - 5 \ 

V l>r6=0 

Since f(S) verifies (14), /(5)e M<5 is still exponentially decaying. The same 
holds for q(5)e~^ s , since q(S) grows polynomially. In turn, this implies that 

?(0=Ee" MS #)/(' + l-5)e M ' +W) 

<5=0 

is the convolution of two exponentially decaying functions, and consequently is 
still exponentially decaying. This completes the proof, with g(r) = X^rMO- D 

5.2 LTQO and frustration freeness 

If £ verifies LTQO and frustration freeness, then local observables are not able 
to distinguish between a global periodic point and a local one. 

Lemma 5.4. Let O a € A a be an observable, po e p\ e Xc A{l) ■ Then 

tx\p OA-PiOA\^\\O A \\K\A\)^)- (17) 
Proof. Frustation-freeness implies that po e Xc A{li) , and thus by LTQO we have 

\\pt-p\<k{\A\)\ {i), 
which implies the desired result. □ 

5.3 Local rapid mixing 

Local rapid mixing is the intermediate step we are going through in order to 
prove theorem 4.6. 
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Definition 5.5 (Local rapid mixing). Take Ac A, and define the contraction 
of T t relative to A as 

n A (T t )= sup sup d tr {T t (p)0 A ,T t T^p)0 A ) 

pZO A eA A 
trp=l ||O a |]=1 

= sup d t r(tv A cT t (p),tT A cT t T (j> (p)). (18) 
trp=l 

We say that L verifies local rapid mixing if, for each A c A, we have that 

V A (T t )$k(\A\)g(t), (19) 
where k(r) grows polynomially in r and g(t) is a fast-decaying function. 

Observation 5.6. It follows from the definition that rj A (T t ) ^ rj B (T t ) whenever 
Ac B. In particular, n A (T t ) s$ n(T t ). 

Note that, in contrast with definitions 4.3 and 4.4, the quantity rj A (T t ) 
depends on the evolution on the whole system, and not just on the subset A. 
Thus local rapid mixing is a very strong condition: the term k(r) appearing 
in equation (19) only depends on the support of A, and so the reduced mixing 
time (i.e. the time it takes for the reduced density matrix on the subset A to 
converge) is independent of system size. On the other hand, the requirement on 
the convergence rate is weaker; it must be faster than any polynomial, but not 
necessary exponential. 

Example 5.7. A simple dissipative system verifying definition 5.5 is the tensor 
product of amplitude damping channels acting (with the same rate) on different 
qubits. Note that, though it might seem a trivial example, there are interesting 
dissipative systems of this form: among others, dissipative preparation of graph 
states [27] can be brought into this form by a non-local unitary rotation (which 
of course does not change the convergence rates). 

Interestingly, if C verifies local rapid mixing, then its fixed points have a 
particular character: they verify a fast-decay of correlations, meaning that the 
correlations between two spatially separated observables decays in distance faster 
than any polynomial. More specifically, if C verifies local rapid mixing with g(t) 
as a decaying function, then the correlations in its fixed points decay as fast as 
g(t). We will suppose for simplicity that T t do not have periodic points beyond 
the fixed points (i.e. = J-£) and that its fixed point is unique. 

Proposition 5.8. Let O a and Ob be operators supported on regions A and B, 
with \\Oa\\ = II Ob || = 1, and let d A B be the distance between them. Define the 
correlation function 

cott(O a ,O b ) = \{O a O b ) - (0 A ) (O b )\ » 

where (O) = tr(/0 oo O) is the expectation value of the operator O acting on the 
fixed point . Then we have that 

corr(CU, Ob) < p(\A\ + \B\)g{d AB ), 

where g(d) is a fast- decaying function and p(-) is growing polynomially. Moreover, 
if C verifies local rapid mixing with an exponential convergence rate, then g(d) 
is also exponentially decaying. 
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Proof. Let us call C = Au B. Fix a product state po (i.e. pq = ® |0)(0|) and call 
p(t) = T t (po). From the definition of local contraction, we have that 

tr [(p(t) - Poo )0 A B ] * r, c (T t ) < fc(lA| + 

where we have used local rapid mixing to bound r\ ' . 
We now use lemma 5.3 to localize the evolution: 

tv[p(t)0 A B ] = tv[p T*(O A O B )] < 

tr [po(T^ s) Y(0 A B )] + \AuB\ ^h(t)f(s), 

where h(t) grows exponentially in time and /(s) is exponentially decaying. If 
s < d AB /2, then T t C(s) is a product of terms acting on disjoint regions, and thus 
we have that T*' C(s) (O a O b ) = T*' A(s) {0 A )T* B{s) (0 B ). 

By applying lemma 5.3 and local rapid mixing one last time, since po is a 
product state we have that 

tr[p a r t ^ s \o A )r t ' B ^\o B )\ 

< (O a ) (O b ) + (\A\ + |S|)-/i(t)/( S ) + [fe(|A|) + fc(|5|)]ff(t). 
Putting everything together, for s ^ d AB /2 we have shown that 

corr(0 A ,O s ) < p(\A\ + \B\)[g(t) + -h(t)f(s)] Vt ^ 0, 

v 

and by choosing a t = 0(d AB ) such that h(t)f(s) is exponentially decaying, we 
can finally bound the correlation as follows: 

corr(0 A , O b ) <p(\A\ + \B\)g'(d AB ), 

where g' is fast-decaying. Note that if g(t) was exponentially decaying, then 
g'(d) is too. □ 

In the following sections, we will prove that the assumptions on the Lind- 
bladians in the stability theorem 4.6 imply local rapid mixing. Together with 
the above result, this implies fast decay of correlations for local Lindbladians 
satisfying these assumptions, as stated in theorem 4.10. 

6 Proof of main result 

6.1 Step la: from truncated to local rapid mixing 

In this section, we want to show that if C verifies frustration-freeness, LTQO, 
and truncated rapid mixing, then it verifies local rapid mixing. This is half of the 
first step needed to prove theorem 4.6; the second half of this first step, involving 
global rapid mixing and uniform families of Lindbladians, will be presented in 
the next section. 

Let us recall a result from [27]. 
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Theorem 6.1 (Contraction for commuting Liouvillians). Let {£j}™ =0 be a set 
of commuting Lindbladians. Define C = Y,j £j a- n d the corresponding evolutions 
T 3 t =e tc i andT t = e tc . Then: 

n{T t )^{Ti). (20) 

3 

In particular, this means that for a local Lindbladian C = £ u C u with finite- 
range interactions composed of commuting terms [£ u ,£ u ] = 0, and such that the 
local terms are uniformly gapped (i.e. rj(T^) s$ ce~ 7 * for some c and 7 independent 
of u and of the global system size) , it holds that 

and thus T t verifies truncated rapid mixing. 

Corollary 6.2. Let £ be a local Lindbladian, with finite-range interactions and 
commuting terms. If the local terms of C are uniformly gapped, then C verifies 
truncated rapid mixing. 

Interestingly, theorem 6.1 also implies that if equation (9) holds for balls, it 
also holds for disjoint unions of balls, since for any disjoint sets A and B the 
respective truncated Lindbladians commute: 

AnB = ^[£ Al £ B ] = 0, 

and consequently r/(T^ uB ) < rj{T A ) + r/(T B ). Thus verifying equation (9) for 
balls is enough. 

We are now ready to prove the key lemma for this step of the proof. 

Lemma 6.3. Let C verify LTQO and frustation-freeness. Then, for all Ac A 
and all i > 0, we have 

V A (T t )4-\A\(e vt -l-vt) g ( S ) +P (\A\)A (s) + n A (T t A(s) ) V,s > 0, (21) 

where g(s) is an exponentially decaying function, p(-) is growing polynomially 
and the constant v is the Lieb-Robinson speed for T t . 

Proof. Let O a be an observable supported on A with ||Oa|| = 1- Fix s ^ 0, and 
call B = A(s). Then, by the triangle inequality 

\\t;(o a )-t;t;(o a )\\ $ \\t;(o a ) - (T t B y(o A )\\ 

+ \\(T t B y(o A )-(T t B T*y(o A )\\ 

+ \\(T t B T B y(0 A )-T;T t *(0 A )\\. 

We can bound the first term in the sum using lemma 5.3 (which used Lieb- 
Robinson bounds to localize the evolution): 

\\T*(O a ) ~ (T B y(0 A )\\ * \A\ -(e vt - 1 - vt)g(s), g exponentially decaying. 

v 

The second term is bounded by the contraction of the localized evolution T B : 
\{T B y{0 A ) - (T t B T B y(0 A )\\ $ V A (T B ). 
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To bound the third term, we rewrite it as follows: 

I (T t B T*) * (0 A ) - T;T*{O a )\\ = sup tr \[T t B T*(p) - T t T^p)]0 A \ . 

p 

Observe that Xc (resp. Xc B ) is closed under the action of T t (resp. T t B ) (and 
the evolution on these spaces is unitary). This means that po '•= T B T B (p) e Xc B 
and pi ■■= T t T^(p) e X c . We can now apply lemma 5.4 (which used LTQO and 
frustration-freeness) to obtain 

tr\(T t B T*(p) -T t T^(p))0 A \ = tr|( Pl - p )O A \ 4 p(\A\) A ( S ) , 

which completes the proof. □ 

Proposition 6.4 (From truncated to local rapid mixing). Let C verify LTQO 
and frustation-freeness. Then, if £ verifies truncated rapid mixing, it also verifies 
local rapid mixing. 

Proof. Fix A c A and s ^ 0. Remember that truncated rapid mixing is a bound 
on the contraction of the local evolution T t : 

f7(Z?)<fc(|4) e -* 

In definition 4.3 this is required to hold only for ball-shaped domains, and we 
know that thanks to theorem 6.1 it extends to disjoint unions of balls. The 
problem we face now, if we want to plug this inequality into the result of 
lemma G.3, is that if A is not connected then A(s) is not of this form when s 
is bigger than the distance between the connected components of A; instead, it 
has a "peanut-like" shape. 

The key observation is that, when two disconnected components of A grow 
enough to "touch" each other, the ball containing both has size polynomial in 
the distance d between them, and thus asymptotically we do not loose anything 
by changing from A(d) to the smallest ball containing A(d). While in general 
?y(T t A ) is not monotone as A increases, Lieb- Robinson bounds and lemma 5.4 
are (in the sense that replacing A(d) by a ball containing it does not affect the 
resulting bounds - it would actually improve them), so it is enough to allow for 
discontinuous steps in the growth of A(s) in order to ensure it is a disjoint union 
of balls. 

Given this, we have that for each s ^ 0: 

V A (T t ) 4 - \A\ (e vt - 1 - vt)g(s) +p(\A\)A (s) + k(\A\)h(s)e-*> , 

v 

where h(s) is bounded by a polynomial in s. 

We want to show that, for each time t 2 0, we can choose s = s(t) in such 
a way that both e vt g o s(t) and e~' 7 h o s(t) have super-polynomial decay in t. 
For each t Z 0, let us choose an s = so(t) such that 

(e vt -l-vt)g(s)^e- vt Vs^s . 

We can choose So(t) such that it grows at most linearly: so(t) = 0(t). Let us 
choose a second Si = Si(t) such that 

h{s)e- lt e"2* Vs^si. 
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We can choose s\(t) such that log(si(i)) = 0(t). 

Since so(t) grows at most linearly in s, while si(t) grows at least exponentially 
in t, by choosing s = So(t) we have that, up to constants which do not depend 
on the system size, s(t) 4 si(t) for all t ^ 0. 

With this choice of s(t), A (s(i)) is still fast-decaying in t, thus 

r, A {T t ) < - \A\e- vt + fc(l4)A„(*(t)) +p(\A\)e~ t i < g(L4|Mi), 

where <p(t) is a fast-decaying function of t and g(-) grows polynomially. This 
completes the proof. □ 



6.2 Step lb: from global to local rapid mixing 

We will now focus on the other set of hypotheses of theorem 4.6, in which the 
truncated rapid mixing hypothesis is replaced by a simpler global rapid mixing 
condition and an assumption on the "self-similarity" of the family of Lindbladians 
being considered. To explain the main idea, consider a translationally-invariant 
model with open boundary conditions. Then the truncated dynamics T^"^' is 
exactly the same as the global dynamics defined on the smaller lattice b x (r). 
This implies that, just by looking at how the global contraction scales with the 
lattice size, we can check the truncated rapid mixing assumption. This scaling 
of the mixing time with the system size is a more natural (and simpler) property 
to consider. 

We would like to generalize this idea to models with arbitrary boundary 
conditions. We first need to pin down some minimal properties that the family of 
Lindbladians should satisfy, for the scaling of their mixing times to be meaningful. 



Definition 6.5. Consider a family of Lindbladians {C l }l^ with finite range or 
exponentially decaying interactions, where each C L is defined on Al c iP such 
that b (L) c A L . We say that the family is uniform if the following conditions 
hold: 

(i) C L are uniformly bounded in strength: 

ll^uMlL , , 

II r-Ll \ II t 1 1 l-»l.c& j,/ x 

sup ||£ tl (r)|| 1 ^ 1 cb = J < oo, sup = f{r), 

L,u,r L.u J 

where f(r) is compactly supported (finite range interactions) or exponen- 
tially decaying. 

(ii) There exists a translationally-invariant 7 Lindbladian M., defined on iP ' , 
such that for all L > 0: 

where dA^ is defined as 

5A i = A i x6 (L). 

7 We restrict to translationally-invariant Lindladians for simplicity, and because this is the 
most important case. The results in fact extend to more general families of Lindbladians, for 
which M in definition 6.5 is not translationally-invariant. 
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In other words, the local terms of the Lindbladians cannot depend on the 
system size, up to boundary conditions (i.e. terms whose support intersects 
8Al), which are allowed to change freely between the different system sizes as 
long as they can be uniformly bounded in norm. This includes, for example, all 
translationally-invariant Lindbladians with periodic, open, or fixed boundary 
conditions. 

Since we are interested in observables whose support is not connected, we will 
use the same convention laid out in the proof of proposition 6.4: with an abuse of 
notation, for a finite A c iP we will denote by A(s) the smallest disjoint union 
of balls containing the set {x\d(x,A) ^ s}. As noted in the previous section, 
|A(s)| still grows polynomially with s, and this notation allows us to define the 
following: 

Definition 6.6. Let A c TL D be a finite set, and let A(s) = \Ji=ib Xi (Li), for 
some Xi e 1j D and some Li ^ 0. We will denote by C A ^ the formal sum 

i=l 

acting on Aj, j . 

If we remove the boundaries of each , then we can identify what is left with 
A(s), by translating each bo(Li) to b Xi (Li), and consequently we can identify 
^A(s) ^' e ' interaction terms of C A ^ which do not touch any boundary) 
with A4a( s ), which in turn coincides with £-a(s) ^ or a ^ ^ ^ diamA(s). 

The following result is the analogue of lemma 6.3. 

Lemma 6.7. Let {£ l }l be a uniform family of Lindbladians verifying LTQO 
and frustation-freeness. Let A c X D be a finite region. Take s ^ and L ^ 
diamyl(s). Then for t ^ 0, we have that 

n A (T^)^^\A\(e^-l-vt)g( S ) + k(\A\)A (s) +V A (T t A{f)s) ) V S ^ 0, (22) 

where g(s) is an exponentially decaying function, fc(-) grows polynomially, the 

T t r ^ A( 3 s} 

constant v is the Lieb-Robinson speed for T t = e , T t is the evolution 
generated by C A ^ S \ and /3 > 1 is a constant such that 

\A((3s)\>2\A(s)\. 

Proof. Let Oa be an observable supported on A with \\Oa\\ = 1- Call St the 
evolution generated by C A ^ S ^ and V t the one generated by = ^a(s)- 

For notational simplicity, we drop the superscript L from T t L . Because we are 
supposing that C L is a uniform family, Vt can be seen as the truncation of either 
T t or S t . And since (3 is such that \A(/3s)\ =s 2 \A(s)\, the 1 ^ 1 norm of V t - V t V<p 
is the same in both cases, because of observation 3.1. 

Following the same argument as in the proof of lemma 6.3, we have that: 

\\T t *(0 A )-T t T;(0 A )\\ 

^ l {e vt _ vt _ 1)g{s) +p{ \a\)a (s) + \\v t *(o A ) - v;v;(o A )\\ . (23) 
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On the other hand, by the triangle inequality we have: 

\\v t *(o A ) - v;v;(o A )\\ $ \\v t \o A ) - s* t (o A )\\ 

+ \\s;s;(o A )-v t *v;(o A )\\ 
+ \\s:(o A )-s;s;(o A )\\. 

Since V* is a truncation of St , we can bound the first term of the r.h.s. by lemma 
5.3, the second term using LTQO, and the third one simply with the contraction 
of S t . 

\\V t *(0 A ) - V t V;{0 A )\\ < -(e vt - vt - l)g(s) + p(\A\)A (s) + r, A (S t ). (24) 

Putting together equations (23) and (24), we get the desired bound. □ 

Finally, we have the analogous result to proposition 6.4. We omit the proof, 
since it follows the same steps as the proof of proposition 6.4, using lemma 6.7 
instead of lemma 6.3. 

Proposition 6.8 (From global to local rapid mixing). 

Let be a uniform family of Linbdladians verifying LTQO and frustation- 

freeness. Then, if it verifies global rapid mixing, 8 it also verifies local rapid 
mixing. 

6.3 LTQO and unique fixed point 

In the hypotheses of theorem 4.6, we did not assume LTQO for the case of 
uniform families of Lindbladians, arguing that it was not needed. We will now 
show that LTQO is already implied by the other assumptions in the theorem. 
For open boundary conditions this is trivial; in this case the truncated dynamics 
are themselves in the uniform family, and the unique fixed point assumption 
implies that there are no states to distinguish in definition 4.2. However, for 
other boundary conditions the result is non-trivial, and follows from the global 
rapid maxing condition. 

Theorem 6.9. Let {C l }l be a uniform family verifying global rapid mixing and 
frustration- freeness, and suppose each C L has a unique fixed point and no other 
periodic points. Then C L verifies LTQO. 

Proof. Fix a finite set A c iP . Let s Z 0, and fix L ^ diam^4(s). Let us denote 
by S t the evolution generated by L A ^ as in definition 6.6, and V t = exp(i£^|^). 

Call poo the unique fixed point of C A< ^ S \ Consider a periodic point p for 
Then by lemma 5.2, for each observable O a supported on A, we have that 

\\V t *{0 A ) - S* t {0 A )\\ ^ \\O a \\ \A\p(s)e-^e vt , 

where p(s) grows polynomially in s. This implies that, by calling p\ = V t ~ 1 (po): 

tr\O A (p -St(p 1 ))\ = tr\O A [Vt(p 1 )-St( Pl )]\$\\O A \\ \A\p{s)e^ s e vt . 

8 If we allow non-translationally-invariant Lindbladians in definition 6.5, we can no longer 
argue that rapid mixing holds for arbitrary balls by translation. Instead, we must require 
global rapid mixing for all Lindbladians obtained by truncating M to an arbitrary finite region 
and adding boundary terms. 
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On the other hand, global rapid mixing implies that 

tr|0 A (S t (pi) - Poo )\ 4 \\0 A \\ n(S t ) 4 \\0 A \\ q(s) \A\ e*, 
with q(s) growing polynomially. Then, by choosing vt = /is/2, we have that 

A (a) = \A\p(a)e-> t 'e vt + q(8)\A\e~T t 
is a fast-decaying function, and 

ti\0 A (po-p«>)\ < \\0 A \\ Ao(s). 

Since pa was an arbitrary periodic point of by the triangle inequality 

(comparing each local fixed point with p^) we obtain LTQO for C A ^\ 

Now let p be a fixed point of ^a(s)- By frustation freeness, we have that 
pA{s) _ tr A ^ s ^ c p i s a fixed point of £^( s ) = an d thus 

\\p A -pi\\^A Q (s). 

□ 

Interestingly, a similar proof also holds when we replace global rapid mixing 
with local rapid mixing. 

Theorem 6.10. Let C verify local rapid mixing, and suppose it has a unique 
fixed point. Then C verifies LTQO. 

Proof. Let poo be the unique global fixed point of C. Consider a region Ac A, 
and consider two fixed points Pq,P\ of Ca(£)- To prove that po and pi arc 
indistinguishable on A, we will prove that both are indistinguishable from 
p^ = tr A c p x , and by triangle inequality 

I A All , II A A\\ II A A I 
\\P0 ~ Pi 111 \\P0 ~ Poo 1)! + \\Pl ~ Poo Id- 

Let Oa be an observable supported on A such that \\Oa\\ = 1- Fix p e Xa(£)- 
We want to find a fast-decaying function Aq(£) such that 

\tv[0 A (p- Poo )]\Kp(\A\)A (£), 

for some polynomially growing function p(-). 

For each t ^ denote T t = exp(t£) and St = exp(t£ A (e))- Call pi = S i ^ 1 (p). 
Then, by lemma 5.3, it holds that 

|tr [0 A (T t ( Pl ) - p)]\ = |tr [0 A (T t ( Pl ) - S t ( Pl ))]\ *c\A\ e vt f(£), 

where f(£) is an exponentially decaying function. 

On the other hand, local rapid mixing implies that 

|tr[0 A (T t ( Pl )- Poo )]| = MOAiTM-T^ipiM^iTt) = k(\A\)g(t), 
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where g(t) is fast-decaying in t. 

Putting these last two bounds together, we have that for all t ^ 0: 

|tr [0 A ( P - Poo )]\ 4 |tr [0 A ( P - T t o S;Hp))]\ + |tr [0 A (T t o S;\p) - Peo )]\ 
$c\A\e vt f(£) + k(\A\)g(t). 

By choosing t = t{£) = 0(£) such that \A\ e vt f(£) is decaying in £, we have proved 
the desired result by defining 

A (£) = c\A\e^f(l) + k(\A\)g(t(£)). 

□ 

6.4 Step 2: from local rapid mixing to stability 

We now prove that local rapid mixing alone implies stability. This is the last step 
in the proof of theorem 4.6, as we already proved in the previous sections that 
the hypotheses of theorem 4.6 imply local rapid mixing. However, the following 
result also stands independently: if a system can be shown to verify local rapid 
mixing by other means, it will also be stable. 

Theorem 6.11. Let £ be a local Lindbladian verifying local rapid mixing, and 
having a unique fixed point p^ such that 

T <l> = T 00 = \p 00 ))((l\. 

Then, using the same notation as in theorem J±.6, for all observables O a supported 
on Ac A we have that 

||Oo(0-Oi(t)Hc(|A|)||0 A ||f, (25) 

for some c not depending on the system size, independent of t, and polynomial 
in \A\. 

Proof. Call T t = e tc and S t = e t{c+E) , and recall that O x = e tc * (0 A ) and 
O2 = e*( c +E '(Oji). Let us write the difference Oq - 0\ using the integral 
representation 

t t 
Q {t)-O x (t)= [ St s E*T:(0 A )ds = Y,Z [ Sl s E u (ryT;(0 A )ds, 

u r 

which implies 

||Oo(*)-Oi(*)K££ [ \\E a (ryo Q (s)\\ds, 

u r q 

where we used the fact that St is a contraction. 

Fix a site u and a positive r, and call S = dist(A, b u (r)). We can split the 
integral at a time t (to be fixed later, depending on S). We bound the first part 
of the integral with Lieb-Robinson bounds: 

r r e vt ° - Vtn - 1 
/ \\E u (ryO (s)\\ ds 4 £f(r) \\0 A \\ \A\ e^ 5 - ^ . 

J v 
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Now pick t =to(S) such that 

e -^Q-,tQ-l^_ ¥ 

V 

We can choose t (S) = f^<5 = 0(<5). 

If t s$ £o(<>), then we have bounded the entire integral, and we are done. 
Otherwise, we treat the second part of the integral as follows: 

t t 
f \\E u (ryO (s)\\ds= J \\E u (ry(O ( S )-T^(O A ))\\ds 

t (5) t (S) 



OO 

£f(r)\\0 A \\ f V A (T s )ds 

to(S) 

OO 

£f(r)\\0 A \\k(\A\) f g(s)ds, 



to(S) 

where we used 

E*(r)T^(0 A ) = E* u (r) \l))({ Poo \ \0 A )) = ( Pao \0 A } E* u {r) |1» = 

together with the local rapid mixing hypothesis. 
Since t (5) is linear in 5, we have that 



g (S) = e + J g(s) ds 



to(<5) 

is still a fast-decaying function. Putting the integral back together, we obtain 
t 

f ||25 u (r)*Oo(*)| ds < £ \\O a \\ fci(|A|)/(r)ff(5). 
o 

Returning to the sum, we have proved that 

\\O (t) -Oi(t)|| * h(\A\)£ \\O a \\ ^^/(r)5(dist(A,6 r (u))). 

u r 

It suffices to show that the sum on the r.h.s. is finite (and independent of system 
size). Let us decompose the sum as follows 

££/(r)5(dist(AA(«))) 

u r 

= E E/(O5(0)+E E (E/(r)flf(d-r)+ E /(r)flT(O)) 

dist(w,A)=0 r d>Odist(M,A)=d \r=0 r=d+l / 

E/(r)3r(r-d) + flT(0) E /(0 . 

r=0 r=d+l / 

where = |{u : dist(u, A) = d}| grows polynomially in d. 
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The first term is clearly bounded, since /(r) is summable. Since / and 7j are 
both fast-decaying functions, their discrete convolution f *'g(d) = Er=o /( r )5 r ( r ~ 
d) is also fast-decaying, and consequently summable against any polynomial. 
The same holds for Y, r >d f( r )- This proves that the second term is also bounded. 
Calling 



c(H) = fc 1 (H) 

we can finally derive the claimed bound 



9(0) \A £ / (r) + £ q(d) If * g(d) + ff(0) £ /(r) ) 

r d>0 V r>d I 



□ 



7 Glauber dynamics 

As an example of a non-trivial dynamics for which we can now prove stability 
using our results, we turn to one of the most studied dynamics in classical 
statistical mechanics: Glauber dynamics, a Markov process that samples thermal 
states of local (classical) Hamiltonians on lattices. Apart from being an interesting 
model in itself, it has important applications in Monte-Carlo Markov chain 
algorithms for numerical many-body physics [33] . Determining whether Glauber 
dynamics is stable against noise or errors is therefore an important - and, as far 
as we are aware, still open - question. 

In this section, we present a natural embedding of Glauber dynamics into the 
Linbdlabian setting, showing how this embedded dynamics inherits properties 
from the classical Markov chain. We will then apply the results of section 4 to 
prove, in the appropriate regime, stability of Glauber dynamics. 

We will consider a lattice spin system over Y = TP or Y = (1i/L1i) D , with 
(classical) configuration space of a single spin a finite set S. For simplicity, we 
will consider the case S = { + 1, -1}. For each A c T, we will denote by Oa the 
space of configurations over A, namely S A . A c will denote the complementary of 
A in T, namely Y \ A. 

Definition 7.1. A finite range, translation-invariant potential { Ja]a<zV is a 
family of real functions indexed by the non empty finite subsets of Y verifying 
the following properties: 

1. j A -.n A ^ R. 

2. For all AcY and all x e Y: 

Ja(o-) = Ja+x(v) "if °~(y + x) = v(y) Vy^A. 

3. There exists a positive r > such that J a = if diam A > r, called range of 
interaction. 

Given a finite range, translation-invariant potential, we can define a Hamil- 
tonian for each finite lattice A c Y and each boundary condition r e fl^a by 

AnA*0 

where a x r is the configuration that agrees with a over A and with t over A c . 
For each such Hamiltonian, we define the Gibbs state state as 

lil(a) = (ZJ,)- 1 e X p(-Hl(a)), 
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where Z\~ is a normalizing constant. 9 The convex hull of the set of Gibbs states 
over A will be denoted by G(A): 

G(A) =conv{/4|Tefi A =}- 

Definition 7.2. The Glauber dynamics for a potential J is the Markov process 
on fl\ with the following generator: 

(Qa/)(<t)= Y,cj(x,a)V x f(a), 

xeA 

where V x f(<r) if defined as f(cr x ) - /(it), and a x is the configuration obtained 
by flipping the spin at position x: 



\a(y) ifxty 
l-cr(x) ifx = y. 



The numbers cj^XjCr) are called transition rates and must verify the following 
assumptions: 

1. Positivity and boundedness: There exist positive constants c rn and Cm such 
that: 

< c m < cj(x, a) < cm < oo Vx, a. 

2. Finite range: cj{x,-) depends only on spin values in b r (x). 

3. Translation invariance: for all k eT, 

cj(x,a') = cj(x + k,a) if a' \y) = cr(y + k) Vy. 

4- Detailed balance: for all x e T and all a 

ex P| " S j a(o-) J cj(x,a) = c,j(x,a x )exp\ - £ Ja(o- x ) J . 

These assumptions are sufficient to ensure that Q\ generates a Markov process 
which has the Gibbs states over A as stationary points. 

Definition 7.3. A quantum embedding of the classical Glauber dynamics for a 
potential J is generated by the following Lindblad operators 

L Xiri = y/cj(x,T)) \r) x ){r]\®l, Vie A,V?7er! Mr) ; (26) 

£x,r,(p) = L x , n pL* xr] - ^{p,Cj(x,T}) \T])(V\} ; 

£aW=EE«,,-1{^}. K =Y J [Y J Cj{x 1 a)\\a){a\- (27) 

xeA r] ^ a \ x / 

plus a dephasing channel acting independently and uniformly on all sites xeA: 
D x , o = Vt|0>(0|, D xA = V7|1)<1| , 2?(p)= S E I?*,*pi?x,i-|A|TP- (28) 

xeA i=0,l 



^Following [37], in our notation we have incorporated the usual inverse temperature param- 
eter directly into the potential J. 
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C\ verifies translational invariance because the transition rates cj do, so this 
family of Lindbladians is uniform. 

Observation 7.4. Take an element of the computational basis, and call 

d(a,(3) the Hamming distance between a and f3. Then it holds that 

V(\a)(f3\) = - 1 d(a,f3) \a)(/3\. 

In other words, I? is a Schur multiplier in the computational basis, with Schur 
matrix given by (-jd(a, /3)) a a. 

On the other hand, we have that for all x: 



cj(x,a) (\a x )((3 x \ - \a){/3\) if a\ K{r) = P\ K ( r) , 
r\ ( c j( x i a ) + c j( x > P)) \ a )(P\ otherwise. 



(29) 

Since d(a x , f3 x ) = d(a, /?), [X>, C XyV ] = for all x e A, and in particular T> and 
C\ commute. 

This quantum dissipative system inherits various properties from its classical 
counterpart. 

Definition 7.5. Let /i be a full-rank positive state. Call 

We say that C is in detailed balance [46] with respect to \i ifT^oC = C*o V ^. 

Proposition 7.6. Let be a Gibbs state over A. Then Ca and T> are in 

detailed balance with respect to jU^. 

Proof. Note that I 1 ^ is a Schur multiplier in the computational basis: 

r MX (|?7i)(?72l) = nKm)* ulto)* \m){m\ . 

From the detailed balance condition for the transition rates cj(x,o~), it follows 
that for all x e A, denoting C x = T,n€a b (r) Ac,?7> 

T MX T M i (|77 1 ){7 72 |) 

= 5 vt,vi \ J ^ X,m \ T {r ll ) ) I 7 ' 1 " 772 ' 2 IWwl 

/ i\i xv./ xi Cj(x,T]-L) + CjyX,r]2) , ,, , 



where 



XX 

V1IV2 



I 1 if Vl\b x (.r) =»72k(r) 

1 otherwise. 



To prove detailed balance for T>, note that Schur multipliers commute, thus 
[2?,r M ] = 0. This, together with the fact that T>* = X>, implies that T> is in 
detailed balance w.r.t. □ 
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The above proposition implies that Gibbs states are stationary states for the 
quantum Glauber dynamics. Let us prove that there are no other fixed points 
apart from the classical ones (i.e. states that are diagonal in the computational 
basis). Clearly, T> has all classical states as stationary points. We just have to 
check C\. 

Proposition 7.7. The set of fixed points of Ca is equal to 5(A), the set of 
Gibbs states over A. 

Proof. Let p be a fixed point of Ca- We want to prove that p is diagonal, i.e. 
that it is of the form 

a 

Consider a non-diagonal element |a)(/3|, and suppose a(x) + /3(x) for some 
x e A. Then, from equation (29), we have that for all y e b x (r), 

£y(l«><0|) = -\(cj(y,a) + cj(y,{3)) \a)(f3\ . 

For y £ b x (r), C y is not supported on x, and thus cannot change the configuration 
there. This implies that the evolution cannot change the configurations over 
the set A(r), where A = {x e A\a(x) + f3(x)}. In turn, this implies that Ca 
commutes with C - Ca (since it acts as a Schur multiplier whose entries depend 
only on the sites in A(r)). Finally, this means that 

h^Hl^WDW, < ||e t£ *(|a>{/3|)| = expf-ti ( £ Cj (x, a) + cj(x, (3) 

V 1 \xeA 

^ exp(-t-c m d(a,P)j -»■ 0. 

Since the off-diagonal elements are killed, p must be of the form Y^aVa l cr )( <J l- 
Writing the equation C\(p) = we obtain 

T l T l c j( x ^ a )p^ \ (jX )(< tX \ ~ EE c ^( x ' cr )^ = °! 

ax ax 

which implies 

Y i cj(x,a x )p cr , = Y.PaCjix.a). 

X X 

The last equation is simply a rewriting of the fact that (p CT ) is a stationary 
distribution for Q\, that is, it is exactly a Gibbs state on A. □ 

Since £a and £a + ^ have the same stationary distributions, even locally, 
all properties that depend just on the structure of the fixed points sets will be 
shared by both: this is the case, for example, of frustration freeness and LTQO, 
which we will prove next. 

Proposition 7.8. C\ (and consequently C^ + T>) is frustration free. 

Proof. By the previous proposition, we have that Xc A = <?(A). We know [33] 
that for Gibbs states it holds that 

AcA^£(A)c£(A), 

but this is exactly the frustration-freeness condition for Ca- Q 
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Let us recall an important definition from [37]. 

Definition 7.9. We say that the Gibbs measures in G(A) satisfy the weak 
mixing condition in V c A if there exist constants C and m such that, for every 
subset A c V, the following holds: 



ytd+v 



-m dist(a;,y) 



T.T'eft 



(30) 



where d^V = {x eV c \ dist(a;, V) s$ r} and fi VA = try N A [iy- 

Proposition 7.10. If Q(A) verifies the weak mixing condition for each V c A, 
then C\ (and consequently C\ + T>) verifies LTQO. 

Proof. Take A c A, I ^ 0, and call V = A(£). The weak mixing condition for V 
implies that there exist constants C and m such that 

sup \^ a -4.a\^C £ e- mdist ^yHCe- me \A\\d + r A(e)\. 

y^d+v 

This is the LTQO condition with A (f) = Ce~ m<! |A| \d+A{l)\. The bound, proven 
for states of the form fi v , can be extended by convexity to all G(V). Let 
VQiVi 6 G(V). By definition, 770 and 771 are convex combination of states of the 
form fiy, thus we can write 

i j i j 

Then we have 



\\Vo,a-Vi,a\\i 



'■■ S VP\\H V ,A- Vv,A\\l 



^ Eft* ~ ^V,a\i 

1 i,j 



□ 



Finally, we want to show that the contraction of the semigroup generated by 
Cj\ + T> can be controlled by the contraction of the classical Glauber dynamics. 
To fix notation, call C ■ A A -* -4a the projector on the computational basis 
diagonal. C is a completely positive, trace preserving map, and it also verifies 
C = lim^oo exp(tT>). Since C\ commutes with V, it also commutes with C. Then 
we can prove the following 

Lemma 7.11. IfT t = exp(t(C A + V)), then 

7 7 (T t )<r7(T t oC) + r 7 (exp(tP)). (31) 
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Proof. Fix an initial state p. Then we can write 

\\T t (p) -TooOOlj * ||T t oC(p) - TMW, + 11^0(1-0(^)11, 

< \\T t o C(p) - ^ oCO^ + lexpCtP) 0(1-0(^)1! 
^r;(T t oC)+7,(exp(i2?)), 

where we have used the fact that £a an d 2? commute, and that the fixed points 
of C\ are invariant under C. □ 

We know, because of theorem 6.1, that 

rKexpC^HIAle-**, (32) 
and this implies the following result. 

Corollary 7.12. If the classical Glauber dynamics verifies global rapid mixing, 
then also the quantum embedded Glauber dynamics generated by C\ + T> does. 

Observation 7.13. Convergence rates of classical Glauber dynamics are a well 
studied subject. It is known that, in some regimes, classical Glauber dynamics 
verifies a Log Sobolev inequality with system-size independent Log Sobolev 
constant (for a review on the subject see [37]). In such situations the classical 
chain has a logarithmic mixing time, and thus verifies global rapid mixing. 

For this class of classical dynamical systems it is possible to apply our main 
result 4.6. In particular, we can arbitrary perturb the transition rates cj(x,o~) 
by some e(x,o~), not necessary preserving detailed balance. If we call £ the 
maximum of |e(x, <t)|, the difference between the perturbed and the original 
evolution of local observables can be bounded by £ times a factor depending on 
the size of the support of the observables taken into account. 

Theorem 7.14. Let the generator of a classical Glauber dynamics, having a 
unique fix point and verifying a Log Sobolev inequality with constant independent 
of system size. Let E be the generator of another classical Markov process of the 
form 

(Ef)(a)=Y.e(x,a)V x f(a). 

Suppose that £ = sup x a \e(x, o~)\ < oo and that e(x,-) has bounded support uni- 
formly in x. Denote by T t the evolution generated by Q\ and by St the evolution 
generated by Qa + E. Then, for each function f supported on A c A, it holds 
that 

\\Tt(f)-S t (f)\\ 00 $c(\A\)\\f\\ oo £, 
for some c(-) independent of system size and polynomially growing. 

Observation 7.15. It is known [35, 36] that the Ising model on 1? or (Z/nZ) has 
a system size independent Log Sobolev constant for high temperatures (when the 
inverse temperature [3 is lower than the critical value /3 C ), or at any temperature 
in presence of an external magnetic field. In this regime the Glauber dynamics 
sampling the Ising model is stable (in the sense of theorem 4.6). 



31 



8 Conclusions and open questions 



We have shown stability of the evolution of local observables under local pertur- 
bations for frustration-free Lindblad evolutions with unique fixed points, under 
the assumption that the mixing time of the system scales logarithmically with 
the system size. In the general case of non-uniform Lindbladians, in which 
no translational invariance is required, additional assumptions of topological 
order and truncated rapid mixing are necessary. We have shown the equivalence 
between LTQO and the uniqueness of the fixed point for such fast-converging 
systems. It should be emphasized that Log Sobolev inequalities provide the 
exact type of convergence time estimates needed for stability. 

The most important open question is what happens outside the rapid mixing 
regime, for example for systems for which the local convergence time depends 
on the global system size. This seems to be the case of state engineering of 
degenerate topological ordered states, such as topologically protected quantum 
codes. 
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Appendix A The non-stable example 

The following example will verify all hypotheses of theorem 4.6, except for 
truncated rapid mixing and belonging to a uniform family, and will be shown to 
be unstable. Interestingly, the system is globally rapid mixing. This example is 
similar to the globally gapped but not locally gapped example in [39]. We will 
show that the characteristics of the dynamics are essentially determined by a 
classical Markov chain embedded into the Lindbladian. For a general review on 
convergence of Markov chains, see [32]. 

Example A.l. Consider a chain of 2N qubits, and define a Hamiltonian H 2N = 
T,1=i Hk,k+i, where Hf-^+i acts only on the qubits (k, k+1) and is defined as 
follows: 

# fe , fe +i = ~U kM1 + |0)(0| fe ® |0)(0| fe+1 + |l)(l| fe ® |l)(l| fc+1 , 

where 

'|l}<l| fc ®|0)(0| fc+1 if k is odd, 



k+1 



if k is even. 

Direct verification shows that H 2N verifies the following properties: 

1. H 2N is local, frustration-free, gapped and verifies LTQO. 

2. |0101 . . . 01) is the unique groundstate. 
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3. 1 1010 ... 10) is the unique lowest-energy excited state, with energy 2/3. 

4. Any other state has energy bigger than 1. 

It is easy to see that this Hamiltonian is unstable under small local pertur- 
bation, beacuse we can perturb it with U = ~Y,k ^Uk,k+i an d the resulting 
Hamiltonian H 2N + U has degenerate groundstate space (we have lowered the 
energy of the first excitation to 0). 

A Lindbladian constructing the groundstate of H 2N is the following: 

2N 2 
fe=l 61S 

where Ck,% are local terms generated by the following corresponding Lindlbad 
operators: if k is odd, then 

i fc ,i=^ +1 |0)<0| fc ®|0)<0| fc+1) 

Lk,2=°* |l)<l| fc ®|l)<l| fc+ i, 
L M =^®^ +1 |1)(1|,®|0)(0|, +1 ; 

if k is even, then 

Lk,i=4 |0>(0| fe ®|0)(0| fe+1 , 
L M = ^ +1 |l)(l| fe ®|l)(l| fe+1 , 
L k ,3 = 0. 

We want to show that C 2N has spectral gap of 1/3. Nonetheless, it fails to verify 
truncated rapid mixing: if we consider truncations of C 2N to subsystems of k 
sites, the corresponding spectral gap is of order O(kfN). 

Consider the elements of the computational basis, which are {|a)(/3|} for each 
a,/3 e {0, 1} 2N . By direct calculation, we can see that the diagonal elements 
{\a)(a\} span a subspace invariant under C 2N . This is the dynamics of a classical 
Markov process that is embedded into our quantum dissipative process. 

Let us analyze the evolution of this invariant subspace first. If we group the 
sites in pairs, forming a chain of N 4-dimensional systems, starting with the pair 
(1, 2), we see that each of the Ck,i can change only the value of one of the pairs. 
This is trivial for k odd; when k is even, Ck.i is supported on two different pairs, 
but acts as a controlled-amplitude damping channel on only one of the two pairs. 
We can then easily write down the transition matrix for the diagonal states: 
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|10) 


|00) 


111) 


|01) 


|io> 
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|00) 
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1 + 1 


111) 










1 + 1 


|01) 


I o 








o ) 



£fc-i, 



Figure 1: The section of transition matrix for C 2N on the pair (k,k + 1) with k 
odd. The blue and the red transitions are present depending on the nearby sites: 
the blue ones if there is a on the right, the red ones if there is a 1 on the left. 
Asterisks in the diagonal are such that the sum of each row is zero. 

Figure 1 shows a local section of the transition graph and its corresponding 
infinitesimal generator matrix. The global graph is given by the cartesian product 
of these sections, and since they are all acyclic (as directed graphs), the global 
graph is acyclic too. To see this, let us define the following 
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The transition matrix for the diagonal states is then given by 



Q 



N 

k=l 
N 

+ Y, ll-..2fe-l 
k=l 

N 



2k-l ® Q ® l2fe+l,...2JV 
Qb® S b ® l2fe+3, 



,2N 



L 2fc+3,...,2Af- 



Therefore Q 2N is upper triangular, the elements on the diagonal are the 
eigenvalues, and consequently are also the eigenvalues of C 2N restricted to the 
diagonal states. We see that the smallest non-zero eigenvalue, corresponding to 
the state |1010 . . . 10), is |. Furthermore, it is easy to see that the diameter of 
the graph of the transitions of Q 2N is N, and this implies that the mixing time 
for C 2N is of order O(logTV) 10 . 



This can be seen from the upper triangular form of Q , noticing that the polynomials 
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Let us now concentrate on off-diagonal elements of the computational basis, 
i.e. matrices of the form |a}(/3] where a * (3 e {0, 1} 2N . Fix one, and call 
A = {k | a(k) * f3(k)}. Then, by direct calculation, we have that 

= -\{p,H kMl } if k or k + 1 e A. 

Given k e A, the only operators that have range over k are and Ck-i, 
but since they only act as Schur multipliers on |a)(/3|, the evolution cannot 
change the value of a and /3 over A. This implies that the terms of C 2N whose 
support intersects A are simply contracting the evolving matrix, while the ones 
supported on A c are acting as the classical Markov chain generated by Q 2N 
(since, by definition of A c , on such sites ]a}(/3] is locally a diagonal element of 
the computational basis). Matrices who share the same set A span an 

invariant subspace for C 2N . 

By permuting the basis in all such invariant subspaces, we can write C 2N fully 
in an upper-triangular form, similarly to what we did before on the "globally 
diagonal" subspace. The eigenvalues of C 2N are then simply given by the 
numbers 

X a , = -^{\aM,H 2N }, Va,/?e{0,l} 2A \ 

and consequently the spectrum of C 2N and the spectrum of H 2N are related as 
follows (where the sum here is intended as a sum of sets) 

2N a(H 2N ) + a(H 2N ) 
a(C )= 2 • 

Since the spectral gap of H 2N is 2/3, the smallest non-zero eigenvalue for C is 
1/3. 

The previous analysis also shows that the fixed points of C 2N are exactly the 
stationary states of Q 2N , and thus C 2N has |01 . . . 01)(01 . . . 01 1 as unique fixed 
point. It is then easy to check that C 2N verifies LTQO and frustation-freeness. 

We have therefore shown that C 2N verifies frustration-freeness, LTQO, and 
is globally gapped. By adding a uniform dephasing noise on all sites T>, by the 
same argument used in lemma 7.11, one can prove that C 2N + V has the same 
fixed points as C 2N , and its convergence rate is given by the classical Markov 
chain Q 2N embedded in the diagonal of C 2N . Since we showed that the global 
mixing time of Q 2N scales as 0(log N), we can conclude that C 2N + T> verifies 
global rapid mixing. 

Nonetheless it fails to verify truncated rapid mixing, since the local interaction 
terms of H 2N are not gapped, and by a similar argument one can prove that 

(f2N , <H 2 -.q k) ) + a { H 2 { l^ 

< J \ J -'(0...,2k)) ~ 2 ' 

Consequently, the spectral gap of £?q 2k ^ is k/3N. 

£2N + 2? is clearly an unstable system, since by removing the terms Ck,3 
(which have strength 0(1/ N)), a new fixed point appears, 1 10 . . . 10)(10 . . . 10] , 
locally orthogonal to the original one. Similarly to example 4.12, the perturbed 
system does verify truncated rapid mixing, but fails to verify LTQO, and is 
obviously unstable. 



appearing in e"2 have degree of at most the diameter of the transition graph. 
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